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ABSTRACT 

' The infrared absorption feature near 2950 cm -1 (3.4 pm), characteristic of dust in the diffuse interstellar 
j medium (ISM), is attributed to C— H stretching vibrations of aliphatic hydrocarbons. We show here that the 
strength of the band does not scale linearly with visual extinction everywhere, but instead increases more 
rapidly for objects near the center of the Galaxy, a behavior that parallels that of the Si— O stretching band 
due to silicate materials in the diffuse ISM, This implies that the grains responsible for the diffuse medium 
aliphatic C— H and silicate Si— O stretching bands are different from those responsible for much of the 
observed visual extinction. It also suggests that the distribution of the carbonaceous component of the diffuse 
ISM is not uniform throughout the Galaxy, but instead may increase in density toward the center of the 
Galaxy, The similar behavior of the C— H and Si— O stretching bands suggests that these two components 
may be coupled, perhaps in the form of silicate-core, organic-mantle grains. Several possible models of the 
distribution of this material are presented and it is demonstrated that the inner parts of the Galaxy has a 
carrier density that is 5 to 35 times higher than in the local ISM. Depending on the model used, the density of 
aliphatic material in the local ISM is found to be about 1 to 2 — CH 3 groups m“ 3 and about 2 to 5 — CH 2 — 
groups m” 3 . These densities are consistent with the strengths of the 2955 and 2925 cm -1 (3.38 and 3.42 pm) 
subfeatures (due to — CH 3 and — CH 2 — groups, respectively) within the overall 2950 cm -1 (3.4 pm) band 
being described by the relations A v j r {295 5 cm ~i> = 270 ± 40 and AJx {2 9 25 cm -i) = 250 + 40 in the local diffuse 
ISM. 

Subject headings: infrared: ISM: lines and bands — ISM: molecules 



1. INTRODUCTION 

The nature of dust in the diffuse interstellar medium (ISM) 
has recently been the subject of a great deal of study. General 
parameters such as particle size and shape have been con- 
strained using the UV-visual portion of the excitation curve 
and polarization measurements. These, in conjunction with 
cosmic abundance constraints and the optical properties of 
relevant materials, have been used to place constraints on the 
dust composition. In general, the dust is believed to consist of 
mixtures of silicate and carbonaceous materials (see Greenberg 
& Hong 1974; Mathis, Rumpl, & Nordsieck 1977; Greenberg 
1978; Greenberg & Chlewicki 1983; Jones, Duley, & Williams 
1987; Mathis & Whiffen 1989). Recent reviews of interstellar 
dust models, and pertinent references can be found in Green- 
berg (1989), Mathis (1989), and Williams (1989). 

The detection of absorption bands near 3300 and 2950 cm -1 
(3.0 and 3.4 pm; due to O— H and C— H stretching vibrations, 
respectively) in the spectra of Galactic center infrared sources 
(Soifer, Russell, & Merrill 1976; Willner et al. 1979; Wick- 
ramasinghe & Allen 1980) led to many follow-up observations 
and analyses because the features along this line of sight were 
thought to be due to dust in the diffuse ISM (Allen & Wick- 
ramasinghe 1981; Willner & Pipher 1982; Jones, Hyland, & 
Allen 1983; Wickramasinghe & Allen 1983; Allamandola 
1984; Butchart et al. 1986; Tielens & Allamandola 1987; 
McFadzean et al. 1989). However, these Galactic center studies 
did not confirm the dust responsible for the features is in the 
diffuse medium, since the features could have potentially been 
produced by materials local to the Galactic center itself. 
Indeed, the work of McFadzean et al. (1989) showed that the 
broad O— H feature near 3300 cm " 1 (3.0 pm) varies in strength 
from source to source in the Galactic center, and they sug- 


gested that it is due, at least in part, to local materials. Without 
additional studies along different lines-of-sight it was difficult 
to ascertain how much of the 3300 and 2950 cm” 1 (3.0 and 3.4 
pm) absorption bands arise from material local to the Galactic 
center and how much (if any) arises from the intervening ISM. 

Until recently, the only other object toward which the 2950 
cm^ 1 (3.4 pm) feature had been detected was VI Cygni 12 
(Adamson, Whittet, & Duley 1990). However, recent studies of 
the 3600-2600 cm” 1 (2.78-3.85 pm) spectral region of objects 
along many different lines of sight suffering from a wide range 
of visual extinctions ( A v ), have clarified this issue considerably 
(Sandford et al. 1991; Pendleton 1994; Pendleton et al. 1994). 
These papers demonstrated that the 3300 cm -1 (3.0 pm) O— H 
stretching band does not correlate with A v along different 
lines-of-sight and therefore is not due primarily to dust in the 
diffuse ISM, confirming the conclusions of McFadzean et al. 
(1989). In contrast, the studies of the 2950 cm” 1 (3.4 pm) C— H 
stretching band by Sandford et al. (1991) and Pendleton et al. 
(1994) show there is remarkable similarity in position and 
profile of this feature over extinctions ranging from as low as 
3.9 to as high as 31 (see Fig, la). In addition, the strength of this 
feature increases with increasing A v , suggesting that it is largely 
due to dust in the diffuse ISM. 

The profile and subpeak positions of the overall C— H 
stretching feature provide specific clues to the nature of the 
carbonaceous material in the diffuse interstellar medium. The 
hydrocarbon feature extends from about 3000 cm -1 (3.33 pm) 
to 2800 cm” 1 (3.57 pm) with subfeatures near 2955, 2925, and 
2870 cm” 1 (3.38, 3.42, and 3.48 pm) (Fig. 1). The positions of 
the first two of these subfeatures are characteristic of the sym- 
metric C— H stretching frequencies of — CH 3 (methyl) and 
— CH 2 — (methylene) groups in saturated aliphatic hydrocar- 
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Fig. 1. — Comparison of the optical depth profile of the C— H stretching 
band due to dust in the diffuse ISM toward GC IRS 6E ( A v — 31 mag ; filled 
circles ) and: (a) VI Cyg 12 (A v = 10 mag; open squares). The VI Cyg 12 data 
have been scaled to have the same optical depth as the GC IRS 6E data at 2925 
cm -1 . ( b ) A C-rich laboratory residue produced by the irradiation of an inter- 
stellar ice analog (solid line) (see Allamandola et al. 1988 for details on the 
production of the laboratory residue), (c) A carbonaceous phase of the Murchi- 
son meteorite ( solid line) [see Pendleton et al. 1994 for further details on the 
comparisons shown in (b) and (c)]. 

bons (molecules with the formula C n H n + 2 )- The relative 
strengths of these two subfeatures indicate that the average 
— CH 2 — /— CH 3 ratio of the interstellar hydrocarbon is 
2.5 ± 0.4. However, the interstellar band profile differs from 
that of pure saturated aliphatic hydrocarbons in that it does 
not contain the two peaks near 2870 and 2850 cm" 1 (3.48 and 
3.51 pm) characteristic of the asymmetric C— H stretching 
vibrations of — CH 3 and — CH 2 — groups, respectively, but 
instead shows only a single subfeature near 2870 cm -1 (3.48 
jnn). Electronegative or other perturbing chemical groups at 
the end of short chain aliphatics can suppress the 2850 cm" 1 
(3.51 ^m) feature without greatly effecting the positions or 
strengths of the other three bands. Thus, comparison of labor- 
atory and astronomical data suggests that the 2950 cm" 1 
(3.4 pm) C— H stretch feature produced by dust in the diffuse 
ISM is largely due to short ( n = 2-4) aliphatic chains 
like -CH 2 -CH 2 -CH 3 and -CH 2 -CH 2 -CH 2 -CH 3 
attached to electronegative or other perturbing chemical 
groups like — O— H, — C = N, aromatics, etc. (see Sandford et 
al. 1991 for a more detailed discussion on the spectral con- 
straints on the diffuse C— H carrier). 


Additional constraints on the nature of this material can be 
provided by spectral comparison with laboratory analogs. In 
Pendleton et al. (1994), we compared the interstellar feature 
with the spectra of refractory carbonaceous residues formed by 
the wanning of UV phot oly zed interstellar ice analogs (see Fig. 
lb ), hydrogenated amorphous carbon (HAC), quenched car- 
bonaceous composite (QCC), and a carbonaceous fraction of 
the primitive meteorite Murchison (see Fig. lc). All these 
materials provided qualitative fits to the interstellar feature, 
but the match provided by the meteoritic material was clearly 
best. This suggests that the perturbing chemical groups 
inferred from the profile of the interstellar feature may consist 
largely of aromatic moieties since meteoritic kerogen is rich in 
aromatics. In this regard, the tentative detection by Pendleton 
et al. (1994) of a weak absorption feature centered near 3030 
cm " 1 (3.3 pm) in the spectra of several Galactic center objects 
may be especially important. This spectral position is charac- 
teristic of the aromatic C“H stretching vibration and suggests 
that aromatic materials may be present in the diffuse dust in 
significant concentrations. 

In this paper we show that, while the optical depth (t ch ) of 
the diffuse medium 2950 cm" 1 (3.4 pm) absorption band 
increases with increasing A v , the observed versus A v rela- 
tion indicates that this material is not uniformly distributed 
throughout the Galaxy. Instead, it appears to be peaked 
toward the Galactic center. 

2. RESULTS AND DISCUSSION 

In the following sections we will demonstrate that (1) the 
carrier of the interstellar C— H stretching feature in the diffuse 
medium is uniformly distributed in the local ISM but is appar- 
ently not uniformly distributed throughout the entire Galaxy 
(§ 2.1), (2) the overall distribution of the carbonaceous carrier 
is very similar to that of silicates in the diffuse ISM (§ 2.2), and 
(3) models that assume an increase in the density of the C“H 
carrier with decreasing galactic radius can explain the observa- 
tions (§ 2.3). 

2.1. The t ch~A v Relation of Dust in the Diffuse ISM 

Since the profile of the diffuse medium C~H stretching 
band is relatively constant from object to object, i.e., the sub- 
features near 2955, 2925, and 2870 cm" 1 (3.38, 3.42, and 3.48 
pm) have very nearly the same relative strengths along all lines 
of sight (Sandford et al. 1991; Pendleton et al. 1994; Fig. la), 
any one of the subfeatures can be used as a tracer of the abun- 
dance of the carrier. We will largely restrict the discussion that 
follows to the 2925 cm" 1 component since it is the strongest, it 
is associated with the more abundant —CH 2 — groups, and it 
does not suffer significantly from spectral confusion due to 
telluric methane. 

Table 1 lists the objects for which we know the strength of 
the interstellar C— H stretching band. Descriptions of these 
objects and their associated visual extinctions can be found in 
Sandford et al. (1991), Pendleton et al. (1994), and references 
therein. With the exception of the objects T629— 5 and OH 
01 —477, which are discussed below, the C— H optical depths 
and visual extinctions given in Table 1 were taken directly 
from Sandford et al. (1991) and Pendleton et al. (1994). All the 
astronomical data used in this paper were obtained at NASA’s 
Infrared Telescope Facility (IRTF). 

OH 01 —477 and T629 — 5 are both M stars whose line-of- 
sight spectra in the 3100-2700 cm" 1 (3.2-3.7 pm) region are 
strongly contaminated by photospheric OH absorption lines 



No. 2, 1995 


GALACTIC ALIPHATIC C-H DISTRIBUTION 


699 


TABLE 1 


Summary of Objects Considered in this Study* 


Object 

Object R.A. 
(1950) 

Object Dec!. 
(1950) 

Extinction 

T b 

''silicate 

t ch 

(2955 cm- 1 ) 

t ch 

(2925 cm -1 ) 

GC IRS 7 

17M2 n 29?3 

— 28°59T3"1 

31 

3.6 + 0.4 

0.19 

0.19 

GC IRS 6E 

17 42 28.9 

-28 59 04 

31 

3.6 ± 0.4 

0.205 

0.21 

GC IRS 3 

17 42 29.1 

-28 5914.9 

31 

3.6 + 0.4 

0.23 

0.24 

T629-5 

18 30 23.0 

-08 44 16.6 

22 



0.07 c 

OH 01-477 

17 46 12.5 

-27 40 49.6 

19 



0.05 c 

AFGL 2179 

18 28 56.8 

-10 01 23.3 

12.8 


0.05 

0.052 

AFGL 2104 

18 13 37.0 

-18 59 49.0 

12.0 


0.037 

0.044 

VI Cygni 12 

20 30 53.4 

+ 41 03 51.6 

10.0 

0.58 + 0.1 

0.043 

0.046 

Ve 2-45 

17 59 00.9 

-23 37 35.2 

6.5 


0.02 

0.02 

BD + 40 4220 

20 30 34.8 

+ 41 08 04 

6.2 


0.027 

0.025 

HD 229059 

20 19 23.0 

+ 37 14 35.0 

5.3 


0.023 

0.028 

AS 320 

18 41 34.0 

-03 5104 

5.2 

0.26 

0.018 

0.02 

HD 194279 

20 21 33.0 

40 35 48.0 

3.9 


0.015 

0.017 


“ Unless otherwise noted, all values are taken from the best available data in Sandford et al. 1991 and Pendleton et al. 
1994. Tch uncertainties are on the order of 20% and are generally dominated by uncertainties in the baselines underlying 
the features, ,4,,-values represent best estimates of the extinction due to the diffuse interstellar material only. All extinc- 
tions are in units of magnitudes. 

b Taken from Roche & Aitken 1985. 

c Values determined after removal of OH photospheric lines. 


which partially obscure the diffuse C— H stretch band toward 
these objects (Sandford et al. 1991). OH 01 —477 is an infrared 
source whose line of sight is very similar to that of the Galactic 
center. This object is coincident with an OH source and 
appears to be a heavily obscured (A v = 19 mag) M4 III star 
(Wickramasinghe & Allen 1980; Jones et al. 1982). T269 — 5 is 
a very late Mira-type, M8 III star reddened by A v « 22 mag 
(Tapia 1981 ; Tapia et al. 1989). 

The C— H band optical depths given in Table 1 for these two 
objects were derived by taking the spectra of Sandford et al. 
(1991) and removing the photospheric OH absorption contri- 
bution by subtracting the spectra of unobscured M stars. For 
unobscured M stars we used ST Ceph and V842 Aql, whose 
spectra were measured by Pendleton et al. (1994). For several 
reasons, this process did not result in optical depths or band 
profiles as precise as those for the other objects listed in Table 
1. First, the data for the obscured and unobscured stars were 
taken during different observing runs and the two sets of data 
do not agree precisely in resolution and detector positioning. 
Second, the absolute and relative strengths of the OH lines in 
M star photospheres are quite sensitive to spectral type, so in 
practice it is difficult to use the spectrum of one M star to 
exactly cancel the OH lines in the spectrum of another M star. 
However, none of the stronger OH lines fall near 2925 cm -1 
(Beer et al. 1972; Sandford et al. 1991; Pendleton et al. 1994) 
and the derived optical depth of the diffuse C— H stretching 
absorption feature at this position is relatively insensitive to 
the interference by photospheric OH, Thus, while it was not 
possible to perfectly remove the photospheric lines and thereby 
determine the profile of the C— H stretching features toward 
T629 — 5 and OH 01 —477, it was possible to use this subtrac- 
tion technique to obtain reasonable estimates of the optical 
depths of the 2925 cm -1 subfeature toward these objects 
(Table 1). 

Figure 2 shows a plot of A v versus t CH(2925 cm -i } for all the 
objects in Table 1 (open and filled circles). The error bars for the 
C— H subfeature optical depths are determined from the larger 
of (1) the statistical error of data points falling in the 2925 
cm' 1 band, or (2) the possible variation of the depth of the 
2925 cm -1 band optical depth when “extreme” flux baselines 


were chosen. In this respect, the error bars are conservative, i.e., 
we show the largest uncertainties. The points near A v = 20 for 
OH 01—477 and T629 — 5 have larger than average error bars 
for the reasons described in the previous paragraph. It is imme- 
diately obvious that the strength of the feature increases with 
increasing A v , a characteristic that would be expected if the 
carrier resides in the diffuse ISM. 

Sandford et al. (1991) demonstrated that the strength of the 
2955 and 2925 cm -1 (3.38 and 3.42 pm) subfeatures observed 
toward the sources they studied could be fit by the relations 
4A<2925 cm -1 ) = 310 ± 90 and A v /r ( 2925 cm -i, = 240 ± 40, the 
largest part of the “uncertainty” arising from the fact that 
these ratios were significantly different in the Galactic center 



Fig. 2. — Visible extinction ( A v ) vs. the optical depth of the C— H and Si— O 
stretching bands of dust in the diffuse ISM. The filled and open circles are the 
optical depths of the 2925 cm -1 (3.42 fim) — CH 2 — band determined from 
low- and high-resolution data, respectively (Sandford et al. 1991 ; Pendleton et 
al. 1994). The open squares are the optical depths of the silicate band near 1000 
cm' 1 (10 (xm) (Roche & Aitken 1985). The silicate optical depths have been 
scaled by a factor of 1/18 to normalize them to the strength of the 2925 cm -1 
— CH 2 — feature at the Galactic Center. (The various points associated with 
the different sources in the Galactic center have been slightly offset from each 
other in A v for clarity). 
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sources studied. These values were subsequently verified by 
Pendleton et al. (1994), who studied additional lines of sight, 
including sources having lower extinctions. They found 
AJx ^ 2955 cm' 1 ) = 270 + 40 and AJx (2 92 s cm- 1 ) = 250 + 40 for 
the local ISM. However, it became clear that the band 
strengths do not increase linearly over the entire A v range, but 
instead rise significantly faster at the higher ^-values associ- 
ated with the Galactic center sources (see Fig. 2 and Pendleton 
et al. 1994). The AJx - values for the Galactic center sources 
alone are AJx {2925 cm -^= 150 ± 15 and AJx l292S cm . l) = 
150 ± 20. Thus, the definition of the diffuse C— H subfea- 
ture strengths by single AJx ratios, i.e., the assumption of a 
uniform distribution of the carrier with A v , may be an over 
simplification. 

At this point it is appropriate to briefly consider whether the 
anomalous AJx ratios for the Galactic center might be due to 
an inaccurate estimate of the visual extinction of these objects. 
For example, if the visual extinction to the Galactic center were 
actually 60 magnitudes instead of 31 magnitudes, the apparent 
anomaly would essentially cease to exist. However, there are a 
number of arguments that suggest that this is not the case. 
Numerous previous studies have addressed the issue of the 
visual extinction towards the Galactic center and the published 
results typically fall in the range of A v = 30-35 magnitudes (see 
Becklin et al. 1978; Henry, DePoy, & Becklin 1984; Sellgren et 
al. 1987; Wade et al. 1987; Rieke, Rieke, & Paul 1989). This 
agreement would seem to be reasonably robust given that 
these values have been obtained from several different types of 
objects requiring different assumed intrinsic colors. In addi- 
tion, using wide field J, //, and K images of the regions sur- 
rounding the Galactic center, Glass, Catchpole, & Whitelock 
(1987) demonstrated that the extinction towards the Galactic 
center is “ patchy ” on large spatial scales due to the distribu- 
tion of foreground interstellar clouds. Examination of the 
maps shows that all three of our Galactic center sources lie at 
the edge of one of these clouds in the H and K maps. Glass et 
al. estimate that the extinction associated with the “ edges ” of 
the clouds seen on their H map correspond to extinctions of 
A v ~ 37 magnitudes. The visibility of the region containing our 
objects on their H maps is therefore in good agreement with 
the published values of A v of 30-35 magnitudes for the Galac- 
tic center sources discussed here. Thus, the anomalous AJx 
ratios associated with the Galactic center appear to be real and 
the assumption of a uniform distribution of the aliphatic C— H 
carrier with A v is not valid. As will be discussed in more detail 
in § 2.3, a more complicated distribution of material seems to 
be implied. 

2.2. A Comparison of the Distribution of Carbonaceous 
Materials and Silicates in the Diffuse ISM 

Included with the x CH -A v relation plotted in Figure 2 ( opened 
and filled circles) is a similar plot of the depth of the 1000 cm -1 
(10 pm) Si— O stretching band of silicates in the diffuse ISM 
(open squares ). The optical depths for the silicate features were 
taken from Roche & Aitken (1985), and have been scaled by a 
factor of 1/18 to match the optical depth of the 2925 cm -1 
band at the Galactic center. The two sets of data are remark- 
ably similar. 

Roche & Aitken (1985) demonstrated that the overall ratio 
of silicate to visual extinction in the Galactic center is about 
twice that of the interstellar medium in the solar neighbor- 
hood. This implies that the grains responsible for the silicate 
band are different from those responsible for the observed 


visual extinction (or at least they are not solely responsible for 
the visual extinction). It also implies that the dust along the 
total line of sight to the Galactic center contains a relatively 
higher proportion of silicates than does the dust in the local 
solar neighborhood (or, alternatively, the Galactic center con- 
tains a lower relative proportion of the material responsible for 
visual extinction). 

The same relative difference between the Galactic center and 
the local ISM is found for the diffuse C— H stretching band 
(Fig. 2). This suggests that the silicate and carbonaceous com- 
ponents of the diffuse interstellar medium may be spatially 
well-correlated and leads to similar conclusions, namely that 
the grains responsible for the diffuse medium C — H aliphatic 
stretching band are different from those responsible for much or 
all of the observed visual extinction and appear to be concen- 
trated toward the Galactic center. In the following section we 
will present three simple models that further quantify these 
points. 

2.3. Models of the Galactic Distribution of Carbonaceous 
Material in the Diffuse ISM 

It is immediately apparent that a simple line cannot fit all 
the data in Figure 2. A good linear fit can be made to all the 
points below A v = 22 mag (solid line in Fig. 2), but the Galactic 
center points at Aj= 31 mag clearly fall well above the 
extrapolation of this line. As mentioned in the previous section, 
the nonlinearity of the x ch versus A v and t siHcate versus A v 
relations indicate that the carriers responsible for the diffuse 
medium C— H and Si— G stretching bands are different from 
those responsible for much of the observed visual extinction. 
Thus, there are two different, simple means of explaining the 
x-A v relations in Figure 2. These are (1) the source of the 
observed visual excitation is distributed throughout the 
Galaxy and the relative abundance of carbonaceous and sili- 
cate materials in the diffuse ISM increases in the inner parts of 
the Galaxy, and (2) the carriers of the C— H and Si— O stretch- 
ing bands are distributed throughout the Galaxy and the rela- 
tive abundance of the source of the visual extinction decreases 
in the inner parts of the Galaxy. Of course, some combination 
of these two, or a more complex distribution, is also a possi- 
bility. In the interest of simplicity, the discussion that follows 
will explore the implications that result if case (1) is true, i.e., 
that the relative abundance of the carriers of the C— H and 
Si — O stretching bands are larger in the Galactic center. The 
reader should keep in mind, however, that nothing in our data 
set preclude the alternative case. We will return briefly to this 
point later. 

Assuming the nonlinearity of the x-A v relations in Figure 2 
results from an increased relative abundance of the C— H and 
Si— O carriers toward the Galactic center, several carrier dis- 
tributions suggest themselves: (1) the distribution of aliphatic 
carbonaceous materials is more-or-less uniform in the “ local ” 
diffuse medium (solar neighborhood) and a smaller volume 
around the Galactic center contains an unusually high concen- 
tration of this material; (2) there is a uniform gradient in the 
density of aliphatic material in the diffuse medium, with the 
density increasing linearly toward the Galactic center; or (3) 
the density distribution of the aliphatic material mimics the 
stellar mass density distribution, i.e., increases exponentially 
toward the Galactic center. These possibilities are explored in 
the following paragraphs using some simple models. 

In all the models we assume the density of — CH 2 — 
absorber responsible for the 2925 cm -1 band (pch 2 ) * s * n un its 
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of — CH 2 ~ groups per cm 2 per unit A v , and we chose M to be 
a variable of extinction running over the range 0 < M < 3 1 
mag. The column density (N*) of any absorbing species, i, can 
be determined from 

N t = J r,{v^dVi/A, * t nmct Av l/2 JA i , 

where J is the integrated area of a band associated with 

the absorber (in cm -1 ), i maxI is the maximum optical depth of 
the band, Av 1/2i is the full width at half-maximum (FWHM) of 
the band (in cm - *), and is the intrinsic strength of that band 
(in cm molecule - *). Thus, for the specific case of the 2925 cm - 1 
(3.42 pm) interstellar — CH 2 — subfeature, we can write 

Ca v 

T CH 2 ~ (^CH 2 /^ V 1/2CH 2 ) J PcH 2 (^f)^M , (1) 

where A CH2 is the integrated absorption strength of the 2925 
cm -1 subfeature due to C— H stretching vibrations in 
— -CH 2 — groups, t CH 2 is the subfeature’s maximum optical 
depth, Av 1/2CH2 is its full width at half-maximum, and p C h 2 (M) 
is the spatial density of the “CH 2 — carrier in groups cm” 2 
mag -1 . 

2.3.1. Model l — A Uniform — CH 2 — Density with a Higher Value 
near the Galactic Center 

In this model it is assumed that the density (p 0 ) of —CH 2 — 
absorber is constant everywhere in the Galaxy except in a 
limited region around the Galactic center, where the material is 
present at a higher density (p GC ). The crossover to the higher 
density domain is assumed to occur at an A v of A 1 where 22 
mag < A t < 31 mag (see Figs. 2 and 3a and Table 2). Thus, we 
have 


pm 


-t 


( 2 ) 


)p 0 for 0 £ M < A Y 
\pgc for A t < M < 31 mag . 
Substituting this form of the density into equation (1) yields 

{ (A C h 2 /Av 1/2 C h 2 )(Po A v) for 0 mag < A v < A t 
( / fcH 2 /AVi/2CH 2 )[(Po ^l) + PgcMv — ^l)] (3) 

for A x < A v < 31 mag , 

which produces an optical depth function (and therefore 


(a) Assumed Gamer Density Function for Model I 








Pi 



0 Aj 31 

Visual Extinction (A v ) 



Fig. 3. — (a) Density distribution of the carbonaceous carrier assumed in 
Model I, i.e., the density (p 0 ) of — CH 2 — absorber is constant everywhere 
except in a limited region around the Galactic center, where the material is 
present at a higher density [p G( f. ( b ) Comparison between the optical depth 
function of the 2925 cm“ 1 band and Model I using the parameter values given 
in Table 2, A x = 25 mag, and p GC = 0.02 — CH 2 — cm -2 mag -1 (see text for 
details). The points have the same meanings as given in Fig. 1. 


column density function) of the form shown by the solid line in 
Figure 3b. 

It is possible to fit the data in Figure 2 with a function of this 
form by making a least squares linear fit to all the i CH data 
below A v = 25 mag (first segment of solid line in Fig. 3b). This 
fit essentially passes through the T629 — 5 point at t CH2 = 0.07 


TABLE 2 


Summary of the Density Models Described in the Text and Their Parameters'* 


Functional Form 

Best- Fit Parameters 

MODEL I: 

^ _ J(Ach 2 I^ v i/ 2 ch 3 )(Po A v ) for 0 ^ A v < A x 

TCH2 1(^ch 2 /Av 1/2CHj )[0? o Aj) 4- p G JA v ~ for A i ^ A v < 31 mag 

N _ JPqA v ) for 0 £ A v < A, 

([(Po^i) + PgJA v ~ AJ] for Ay ^ A v < 31 mag 

p 0 = 0.0032 (Av l/2CH2 /A CH2 ) mag 1 
Pgc = 0.014—0.104 (Av 1/2CHj //I CH2 ) mag 1 

22 mag ^ A x < 31 mag 

MODEL II: 


T CHj “ ( / tcH 2 /^ V l/2CH 2 )fPo A v + ^CA 2 ] 

Po = 0.0012 (Avj^ 2CH2 / 2 4 CHj ) mag 1 
C — 0.00032 (Av j / 2CH2 /A CH2 ) mag 2 

MODEL III: 


t ch 2 — McH 2 /^ v i/ 2 cH 2 )f — C 2 + (p 0 — C 2 )A v + \C x Al 4- C 2 eAJ 
N - j>C 2 4- {p 0 - C 2 )A v 4- \CyA\ 4 C 2 e A '] 

Po = 0.0027 {Avy/jcm/AcM^ mag 1 
C x = 7.16 x 10 5 (Avj mag 2 
C 2 = 2.52 x 10 15 (Av 1/2 ch 2 /^ch 2 ) ma 8 1 


a For the diffuse medium 2925 cm -1 C— H stretching subfeature due to — CH 2 — groups, use A ^ = 8 x 10“ 18 cm/CH 2 
group and Av t/2CHj = 20 cm" 1 (Sandford et al. 1991). All extinctions are in units of magnitudes. 
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TABLE 3 


Summary of the Calculated Local and Galactic Center Densities of 
Aliphatic — CH 2 — in the Diffuse Interstellar Medium* 


Local — CH 2 — Density 
(Po) 

Galactic Center — CH 2 — 
Density (pcj 


(Pgc/Po) 

Model I: 




0.0032 (Av j uchJ^ch) ma S 

0.014 (A v l f 2CH jA CH ^ mag 


4.4 


(A x = 22 mag) 

(-4, 

, = 22 mag) 


to 


to 


0.104 (^v l/2CH JA CH2 ) mag 


32.5 


(A, = 30 mag) 

Mi 

i = 30 mag) 

8.0 x 10 15 CH 2 cm' 2 mag' 1 

3.5 x 10 16 CH 2 cm' 2 mag' 1 


4.4 


{A v = 22 mag) 

Mi 

l = 22 mag) 


to 


to 


2.6 x 10 17 CH 2 cm' 2 mag" 1 


32.5 


{A x - 30 mag) 

Mi 

l = 30 mag) 

5.2 x 10' 6 CH 2 cm" 3 

2.3 x 10' 5 CH 2 cm' 2 


4.4 

= 5.2 CH 2 m" 3 

{A x = 22 mag) 

M, 

[ = 22 mag) 


to 


to 


1.7 x 10" 4 CH 2 cm' 3 


32.5 


(A x — 30 mag) 

Mi 

t = 30 mag) 


= 23-170 CH 2 m" 3 



Model II: 




0.0012 (Av, / 2CH2 /-4 C h j ) mag 

0.01 1 (^v i/2Cli JA CHj ) mag 


9.2 

3.0 x 10 15 CH 2 cm' 2 mag -1 

2.8 x 10 16 CH 2 cm' 2 mag" 1 


9.2 

2.0 x 10' 6 CH 2 cm -3 

1.8 x 10' 4 CH 2 cm" 3 


9.2 

= 2.0 CH 2 m 3 

= 18 CH 2 m" 3 



Model III: 




0.0027 (AvjjjchjMcri) ma S 

0.078 (Av 1/2C>i2 //4 CH2 ) mag 


28.6 

6.9 x 10 15 CH 2 cm' 2 mag' 1 

2.0 x 10 17 CH 2 cm" 2 mag' 1 


28.6 

4.4 x 10“ 6 CH 2 cm" 3 

1.3 x 10" 4 CH 2 cm' 3 


28.6 

= 4.4 CH 2 cm' 3 

= 130 CH 2 cm' 3 




* For the diffuse medium 2925 cm -1 C— H stretching subfeature due to — CH 2 — 
groups, use A clil = 8 x 10' 18 cm/CH 2 group and Av 1/2CH2 = 20 cm -1 (Sandford et al. 
1991). All extinctions are in units of magnitudes. 


and A v = 22 mag. Substituting these values into the top half of 
equation (2) yields p 0 = 0.0032 (Av 1/2CH2 /A CH2 ) ma 8 For the 
2925 cm -1 — CH 2 ~ band, Av 1/2CH2 « 20 cm -1 and A CH2 » 8 
x 10" 18 cm/CH 2 group (Sandford et al. 1991). Assuming an 
average extinction-to-distance ratio of 2 magnitudes per kpe 
(Allen 1976), this corresponds to a local density of ~5 
— CH 2 — groups m -3 (Table 3) and is consistent with a local 
value of A v /T i2925cm -i) = 250. The local diffuse medium 
density of H is ~10 5 H atoms m' 3 (Allen 1976). Since the 
— CH 2 — /“CH 3 ratio of the diffuse medium carrier is 2.0-2.5 
(Sandford et al 1991) this corresponds to a local density of ~2 
to 3 — CH 3 groups m -3 . 

Using equation (3) and the derived value for p 0 , one can 
place limits on the density of aliphatic material implied for the 
Galactic center, where t ch 2 — 0*2 and A v = 31 mag. In this case 
we can only derive limits since we do not know the appropriate 
value of A t . However, A x must lie above A v = 22 mag and 
below A v = 31 mag (Fig. 2). Substituting the appropriate 
values into the lower half of equation (2) yields a lower density 
limit (corresponding to A x « 22 mag) of Pgc = 0.014 
(Av 1/2C h 2 Mch 2 ) mag -1 an and upper limit (corresponding to 
A t * 30 mag) of = 0.104 (Av 1/2CH2 //1 CH2 ) mag" 1 . Using the 
same values of Av 1/2CH2 , A c Hl , and extinction-to-distance ratio 
given above, this corresponds to a range in Galactic center 
densities of about 25 to 170 “CH 2 — groups m' 3 (Table 3). 
Thus, Model 1 requires the density of diffuse medium aliphatic 
materials (— CH 2 “ and — CH 3 groups) to be a factor of ~5 to 
35 times higher in the Galactic center than it is locally, depend- 
ing on the value of A t . 


2.3.2. Model II — A — CH 2 “ Density that Increases 
Linearly toward the Galactic Center 

In Model II it is assumed that the local density of — CH 2 — 
absorber is p 0 and that the density increases linearly as one 
approaches the Galactic center (Fig. 4 a). Thus, we have 

p(M) = p 0 + CM for 0 < M < 31 mag , (4) 

where C is a constant. Substituting this form of the density into 
equation (1) yields 

t ch 2 = (^cH 2 /Av 1/2C H 2 )[p 0 A v + iCA 3 ] , (5) 

which produces an optical depth function (and therefore 
column density function) of the form shown as a solid line in 
Figure 4b. 

The best-fit curve to the data points in Figure 2 having the 
functional form of equation (5) and constrained to pass 
through the origin is = 0.0012A V + 0.00016A 2 (solid line 
in Fig. 46, Table 2). Equating the first-order terms yields a 
local space density of — CH 2 — groups of p 0 = 0.0012 
(Av 1/2CH2 /A CH2 ) mag -1 . This is a factor of —3 lower than the 
local density derived using Model I. Using the same values of 
A v i/ 2 ch 2 > ^ch 2 > an d extinction-to-distance ratio used in Model 
I, this corresponds to a local density of ~2 — CH 2 — groups 
m~ 3 (Table 3) and about 1 —CH 3 group m -3 . 

Equating the second order terms of equation (5) and the best 
fit to the data yields C = 0.00032 (Av 1/2CH2 /Ach 2 ) mag -2 . Sub- 
stituting the appropriate values for p 0 and C into equation (4) 
and using A v = 31 mag yields a calculated Galactic center 
density of p ^ = 0.011 (Av 1/2CH JA CH2 ) mag -1 . This is a factor 
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Fig. 4. — (a) Density distribution of the carbonaceous carrier assumed in 
Model II, i.e., the local density of — CH 2 ~” absorber is p 0 and the density 
increases linearly toward the Galactic center. (6) Comparison between the 
optical depth function of the 2925 cm " 1 band and Model II using the param- 
eter values given in Table 2 (see text for details). The points have the same 
meanings as given in Fig. 1. 

of 1.3 to 9.5 times smaller than that derived using Model I. 
Using the usual values of Av 1/2CH2 , A CH2 ,and extinction-to- 
distance ratio, this corresponds to a Galactic center density of 
^20 —CH 2 ~ groups m” 3 (Table 3) and ~ 10 — CH 3 groups 
m“ 3 . Thus, Model II requires the density of diffuse medium 
aliphatic materials to be a factor of — 10 times higher in the 
Galactic center than it is in the local diffuse ISM. 

2.3.3. Model III — A “CH 2 ~ Density that Increases 
Exponentially toward the Galactic Center 

In Model III it is assumed that the density of — CH 2 — 
increases exponentially as one approaches the Galactic center 
and has a local density of p 0 (Fig. 5a), i.e., 

p{M) = (p 0 — C 2 ) + CjAf + C 2 e M for 0 < M < 31 mag , 

( 6 ) 

where C t and C 2 are constants. This distribution mimics the 
form of the stellar density function in the galactic plane (see 
Allen 1976, p. 285, and references therein). Substituting the 
above form of the density into equation (1) yields 

T CH 2 = (^CH 2 /^ V 1/2CH 2 )[” ^2 + (Po — C 2 )A v 

+ + C 2 e Ao ] , (7) 

which produces an optical depth function (and therefore 
column density function) of the form shown as a solid line in 
Figure 5b. 

The best-fit curve to the data points in Figure 2 having the 
functional form of equation (7) is t CH 2 = —2.52 x 10“ 15 


+ 2.74 x 1( r 3 A v 4- 3.58 x 10 ~ 5 A 2 V + 2.52 x 10” 1 V* (solid 
line in Fig. 5b, Table 2). Equating the zeroth- and second-order 
terms of this expression with the same terms in equation (7) 
yields C x = 7.16 x 10“ 5 (Av l/2CH mag” 2 and C 2 = 2.52 
x 10” 15 (A Vi /2C H 2 /AcH 2 ) mag” . Equating the first order 
terms and using the derived value of C 2 , we find a local 
— CH 2 — group density ofp 0 = 0.0027 (Av 1/2 ch 2 Mch 2 ) mag” 1 . 
This is similar to the local density derived using Model I and 
~2 times higher than in Model II. Using the same values of 
Av iy2 cH 2 > ^ch 2 > and extinction-to-distance ratio used in Model 
I, this corresponds to a local density of ~4-5 — CH 2 — groups 
m ” 3 (Table 3) and about 2 — CH 3 group m ” 3 . 

Substituting the appropriate values for p 0 , C u and C 2 into 
equation (6) and using A v = 3\ mag yields a calculated Galac- 
tic center density of p GC = 0.078 (A v 1 / 2 ch 2 M CH2 ) mag -1 . This is 
a factor of 7 times larger than that derived using Model II and 
falls in the upper part of the range derived using Model I. 
Using the usual values of Av 1/2CH2 , ^4 C h 2 * anc * extinction-to- 
distance ratio, this corresponds to a Galactic center density of 
— 130 — CH 2 — * groups m” 3 (Table 3) and ~ 70 — CH 3 groups 
m” 3 . Thus, Model III requires the density of diffuse medium 
aliphatic materials to be a factor of ~30 times higher in the 
Galactic center than it is in the local diffuse ISM. 

2.3.4. Implications of the Model — CH 2 — Density Distributions 

The parameters of Models I, II, and III, along with their 
predicted local and Galactic center diffuse medium — CH 2 — 
densities, are summarized in Tables 2 and 3. All three models 



Fig. 5. — (a) Density distribution of the carbonaceous carrier assumed in 
Model III, i.e., the local density of — *CH 2 — absorber is p 0 and the density 
increases exponentially as one approaches the Galactic center, (b) Comparison 
between the optical depth function of the 2925 cm ” 1 band and Model III using 
the parameter values given in Table 2 (see text for details). The points have the 
same meanings as given in Fig. 1. 
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indicate that the density of the diffuse medium aliphatic carrier 
is higher in the inner Galaxy than it is in the solar neighbor- 
hood. The extent to which the density in the inner Galaxy is 
higher is model dependent, but it appears that it must be 
higher by at least a factor of 5. Galactic center density enhance- 
ments by factors as large as 35 are required if the local density 
is typical of most of the Galaxy. 

It is presently difficult to choose between these carrier dis- 
tributions because of the paucity of good data at extinctions 
between 20 and 30 magnitudes, the range that discriminates 
most strongly between the models. Models I and III seem to 
produce better fits to the data than Model II, but this state- 
ment is based largely on the better fit to the two least certain 
points near A v = 20 mag. Clearly, additional data in the 20 to 
30 magnitude extinction range would be useful for further con- 
straining the spatial distribution of this carbonaceous carrier 
in the diffuse ISM. 

At this point it is appropriate to remind the reader of several 
points. First, while the previous discussion was addressed spe- 
cifically to the distribution of the carbonaceous material in the 
diffuse ISM, the close match shown in Figure 2 to the diffuse 
silicate distribution determined by Roche & Aitken (1985) indi- 
cates that the same general conclusions apply to this material 
as well, i.e., silicates in the diffuse medium would be similarly 
concentrated in the inner portions of the Galaxy and by similar 
factors. Second, we remind the reader that the above dis- 
cussion was based on the assumption that the nonlinearity of 
the C— H and Si— O stretching band z — A v relationships is 
due to relative abundance enhancements of their carriers with 
respect to the distribution of the source(s) of visual extinction 
toward the center of the Galaxy. Presumably, one could also 
explain Figure 2 by assuming that the C — H and Si — O band 
carriers are distributed throughout the Galaxy and there is an 
under-abundance of nonaliphatic, nonsilicate, extinction- 
producing material in the central regions of the Galaxy. In this 
case, the implied abundance gradients would have similar 
magnitudes to those discussed above, but they would be 
reversed in sense and would apply to the material responsible 
for the visual extinction, not the C — H and Si— O carrier 
materials. 

2.4. Potential Causes of the Nonuniform Distribution of the 
Diffuse ISM Aliphatic C — H Carrier 

The origin of the nonuniformity in the distribution of the 
diffuse medium carriers of the aliphatic C — H and silicate 
Si^O stretching bands with respect to visual extinction is not 
clear. If the actual distribution of the carrier is similar to the 
step function used in Model I, the discontinuity at A v — A i 
(Fig. 3a) might be associated with the transition from the ISM 
in the galactic disk to the central bulge. Alternatively, it might 
be associated with some more localized structure associated 
with the center of the Galaxy, for example, the circumnuclear 
molecular disk known to surround the Galactic center (see 
Kundt 1990 for a review). 

The nonuniformity of the C— H and Si— O carriers may be 
the result of metallicity gradients within the Galaxy. Observa- 
tions of H n regions and supemovae remnants (SNR) in other 
spiral galaxies show that an increase in the O/H, N/H, and S/H 
ratios of about 3 to 5 is not uncommon between the outskirts 
and centers of such galaxies (see Dopita, D’Odorico, & Benve- 
nuti 1980; Blair, Kirshner, & Chevalier 1982; Blair & Kirchner 
1985). Observations of H ii regions and SNR within our own 
Galaxy indicate similar gradients for O/H, N/H, and C/H (see 


Binette et al. 1982; Shaver et al. 1983; Fesen, Blair, & Kirshner 
1985; Meyer 1989). Such gradients appear to be a natural 
result of galactic chemical evolution (see Matteucci & Francois 
1989; Matteucci 1991 ; Rana 1991, and references therein). 

Thus, the increased spatial density of the carbonaceous 
carrier of the aliphatic C— H band toward the Galactic center 
might be the result of an increased relative abundance of avail- 
able carbon. If this were the case, then one would expect the 
behavior of the z CH -A v relation to mimic, to some extent, 
metallicity. In this regards, it is interesting to note that recent 
studies of the abundance gradients of N, O, Ne, and S in our 
Galaxy not only confirm that metallicity increases toward the 
Galactic center, but suggest the possibility that the increases 
may occur in discrete steps rather than continuously (Simpson 
et al. 1995). Such a noncontinuous distribution is consistent 
with the fit provided to our z CH ~A v relation by Model I. 

Unfortunately, this explanation produces a new dilemma, 
namely, it requires that the materials responsible for the visual 
extinction do not track metallicity, or at least not in the same 
way as do the carriers of the C— H and Si— O bands. Given 
that the material responsible for the visual extinction must also 
contain an abundance of elements heavier than H and He, this 
comes as a surprise. Perhaps both types of material track 
metallicity, but the carriers of the C— H and Si— O bands 
increase with metallicity at a much greater rate. 

Finally, in their studies of the Si— O stretching feature of 
silicates in the diffuse ISM, Roche & Aitken (1985) assumed 
that the nonlinearity of their i Si0 -A p relation was due to the 
depletion of carbon stars in the central regions of the Galaxy, 
noting that models like that of Mathis et al. (1977) predict that 
most of the visible extinction comes from grains in carbon-rich 
environments. This type of scenario was mentioned briefly 
earlier and is the compliment of the models discussed here, i.e., 
it corresponds to a model in which the carriers of the C— H 
and Si— O stretching bands are more-or-less uniformly distrib- 
uted in the Galaxy and the relative abundance of the source of 
the extra visual extinction decreases in the inner parts of the 
Galaxy. As mentioned previously, there is nothing in the 
current z-A v data that precludes this possibility. However , this 
interpretation is only consistent with our findings if the carbon- 
aceous material produced by carbon stars is different from the 
carbonaceous material responsible for the 2950 cm -1 band, i.e. 
the C— H band carrier is not made by carbon stars and the 
material made by these stars produces significant visual extinc- 
tion. We note that this possibility is consistent with the inter- 
pretation that the carrier of the diffuse ISM C— H stretching 
band is due to refractory organics produced in dense clouds by 
the irradiation of ice mantles (Sandford et al. 1991), not in 
carbon stars, and offers a natural explanation of the similarities 
in the T S io~^i; and t ch - ^ relations since the silicate grains may 
constitute a large fraction of the core grains on which the ices 
are condensed. In this regard, our findings are consistent with 
the kinds of core-mantle grains suggested by Greenberg and 
discussed in Greenberg (1978) and Hong & Greenberg (1980). 

3. CONCLUSIONS 

The infrared feature in absorption near 2950 cm -1 (3.4 pm) 
in the spectra of many stars suffering from interstellar extinc- 
tion is attributed to aliphatic C— H stretching vibrations in 
carbonaceous dust in the diffuse interstellar medium (ISM). 
The strength of the diffuse C— H stretching band (tch) does not 
scale linearly with visual extinction {A p ) throughout the 
Galaxy, but instead increases more rapidly in the vicinity of the 
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Galactic center. This implies that the grains responsible for the 
diffuse medium C— H aliphatic stretching band are different 
from those responsible for much of the observed visual extinc- 
tion. This also suggests that the distribution of the carbon- 
aceous component of the diffuse ISM is not uniform 
throughout the Galaxy. 

Several simple models of the distribution of the material 
were presented. It was demonstrated that, while the inferred 
increase in carrier density towards the Galactic center is model 
dependent, the inner part of the Galaxy has a relative carrier 
density that is at least 5 times higher than the local ISM and 
which could be as much as 35 times higher. Depending on the 
model used, the local space density of — CH 2 — groups in the 
carrier can range from about 2 to 5 — CH 2 — groups m -3 . The 
density of — CH 3 groups in this material are a factor of 2.0 to 
2.5 times smaller. These densities are consistent with the 
strengths of the 2955 and 2925 cm -1 subfeatures (due to 
— CH 3 and — CH 2 — groups, respectively) being described by 
the relations AJ \ 295$cm ~x ) = 270 ± 40 and Ajx { 29 25 cm-i) = 
250 ± 40 in the local diffuse ISM. Additional high quality mea- 
surements, especially of objects suffering extinctions in the 
A v = 20-30 mag range, will be needed to better constrain the 
carrier density distribution. 
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The relationship between the strength of the 2950 cm -1 (3.4 
fim) diffuse band (t ch ) and visual extinction ( A v ) closely 
matches that of silicate materials in the diffuse ISM, suggesting 
that the silicate and organic components of the diffuse inter- 
stellar medium are spatially well-correlated. In this regard, our 
findings are consistent with the idea that the silicates and the 
carbonaceous material responsible for the 2950 cm -1 band are 
in the form of core-mantle grains. The close correlation of the 
strengths of the silicate Si— O and aliphatic C— H stretching 
bands and their simultaneous lack of linear correlation with A v 
indicates that a substantial portion of the visual extinction in 
the diffuse interstellar medium is produced by an additional, 
spatially uncorrelated component or components. 
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